






Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1960 


Determination of the coefficient of eddy 
viscosity as a function of height and time: and 
separately as a function of stability and 
vertical shear of the wind 


McBride, Jeremiah R. 


Monterey, California: U.S. Naval Postgraduate School 
http://ndl.handle.net/10945/11971 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 


| Calhoun is the Naval Postgraduate School's public access digital repository for 
th D U DLEY research mate = and institutional publications c reated by she NPS community. 
«iit | : 3 Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 
IN | KN Ox appointed — and published -- scholarly author. 
: ] LIBRARY Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


NPS ARCHIVE 
1960 





MCBRIDF, J. 


DETERMINATION OF THE GOEFFIGIENT OF 
EDDY VISCOSITY AS A FUNCTION OF HEIGHT 
AND TIME; AND. SEPARATELY AS A FUNCTION 

OF STABILITY AND VERTICAL SHEAR OF THE WIND 


JEREMIAH R. MeBRIDE 


Library 
U. S. Naval Postgraduate School 


Monterey, California 


PUBLEY KNOX LIRRARY 


NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CA 93943.8704 

















DETERMINATION OF THE COEFFICIENT 
OF EDDY VISCOSITY AS A FUNCTION 

OF HEIGHT AND TIME; AND SEPARATELY 
AS A FUNCTION OF STABILITY AND 


VERTICAL SHEAR OF THE WIND 


RAISE 


Jeremiah R. McBride 





DETERMINATION OF THE COEFFICIENT 
OF EDDY VISCOSITY AS A FUNCTION 

OF HEIGHT AND TIME; AND SEPARATELY 
AS A FUNCTION OF STABILITY AND 


VERTICAL SHEAR OF THE WIND 


by 


Jeremiah R. McBride 
yy), 


Lieutenant, United States Navy 


Submitted in partial fulfillment of 
the requirements for the degree of 


MASTER OF SCIENCE 
IN 
METEOROLOGY 
United States Naval Postgraduate School 
Monterey, California 


1960 





DETERMINATION OF THE COEFFICIENT 
OF EDDY VISCOSITY AS A FUNCTION 

OF HEIGHT AND TIME; AND SEPARATELY 
AS A FUNCTION OF STABILITY AND 


VERTICAL SHEAR OF THE WIND 


by 


Jeremiah R. McBride 


This work is accepted as fulfilling 
the thesis requirements for the degree of 
MASTER OF SCIENCE 
IN 
METEOROLOGY 
from the 


United States Naval Postgraduate School 





ee eee.) 





ABSTRACT 


A harmonic analysis of the observed diurnal variation of the 
coefficient of eddy viscosity is conducted at selected levels followed 
by a mathematical analysis of the amplitude coefficients in order to 
determine a function of height and time needed to represent this 
variation, It is found by this approach that a lengthy series of terms, 
consisting of exponential and polynomial functions of height and 
trigonometric functions of time and height, would be needed to represent 
the observed variation. 

A two layer empirical model is developed from the observed data 
which gives the coefficient of eddy viscosity as a function of stability 
and height in the lower layer, and as a function of vertical shear of the 
horizontal wind and height in the upper layer. 

The author wishes to express his appreciation for the assistance and 
encouragement of Professors G. J. Haltiner and A. B. Mewborn of the 


U. S. Naval Postgraduate School in this investigation. 
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ey Imeroduccson, 

In-an, atcemp. to @ccount> for the observed wind) scructure ink che 
friction layer, meteorologists have devoted considerable effort to 
determine che variation of the coefficient of eddy viscosity with 
respect to time and space. 

The layer of frictional influence has been subdivided into the 
laminar, surface, and spiral layers each possessing certain 
characteristic properties. The coefficient of eddy viscosity, M ; 
has been shown to vary linearly with height in the neutral surface 
layer and as a power function of height in the non-neutral surface 
layer. Using wind data obtained from double theodolite observations 
and assuming a steady state, lildner L6] solved the horizontal 
equation of motion for M utilizing finite differences for the 
derivatives of velocity. He showed that AX increased up to a fairly 
low elevation in the spiral layer (250 meters) and then decreased to 
a small residual value near the cradient level. 

The Taylor Las solution for the wind in che Gag. cr 2535 ec 
KI constant, along with a number of other broad ecu ens, and 
yields a feneral representation of the observed vertical wind variation. 
More recently, Buajitti and Blackadar [ 2| and Haltiner | 4 | have 
sought the distribution of eddy viscosity as a function of time and 
height in order to account for the diurnal variations of the wind 


structure in the friction layer. Haltiner La | assumes the 


coefficient of eddy viscosity to have the form 


Ke HM = g(t) k(z) 


a 





wnere 


e(t) = 1 4b (Sinw +4 c Sin 2W t) 


Cc 


and k(z) takes on several different exponential forms. 

It is the purpose of this investigation to determine the degree 
of complexity of the function oe 4 (z, +) needed to yield observed 
values ogy fa ; and also to determine, if possible, the coefficient 
of eddy viscosity as a function of stability and vertical shear of the 


wind r) 





2% The Coefficient of Eddy Viscosity as a Function of Height and Time. 
The horizontal equations of motion for turbulent flow where the eddy 
stresses are assumed to be a function of height and time only and the 


molecular or viscous stress is considered negligible are 


° ‘\ 3 pret ries 

ee ¢(v- “Va) + 5 

oe ay oT; 

x + Vys -£ (4g) 5 St i 


Here ley and Tzy are the components of eddy stress; 

Ua and a, the components of the geostrophic wind; V » Ve Meaw 

wind velocity vector of whichU andV are the horizontal components; 
ra » the density; and a ; the coriolis parameter. 


The eddy stresses are assumed to be related to the mean wind shear 


by the expressions 


aeons 
ee oe Jey=h < (2) 


The geostrophic departure method employs observations of the mean 


velocity and geostrovhic wind to evaluate all the terms of eq (1), so 


pemeo yield the distribution of oe ont OTe «ais. mernod 
OZ 
fives the distribution of stress and eddy viscosity provided the wind 


uu 
components have a maximum at some levels h, and hy so that at Bx? a 


Ogg 
and at hy » S72 = 


by mulciplying ».é¢q (1) by ? and integrating from h, to an 
arbitrary height 2, the distribution of l7x is obtained. Similarly 
the distribution of |zy is found. Then ea.(2) leads to two independent 


estimates of the distribution of eddy viscosity. 





f ~ s Zz. 
ned I abe J4 = me 
K : 


Assuming no discrepancies in the procedure or dates Ae =u Oe 


Blackadar [1] used the geostrophic departure method on a 
composite wind distribution obtained by averaging the winds at two 
stations (Oklahoma City, Oklahoma and Wichita, Kansas) over 24 
selected, not necessarily consecutive, one-day periods. Figures 1 
and 2, reproduced from Blackadar's[_1} data, show the eastward 
and northward components of the average winds during the 24 day 
period, Although aeaat care was taken in determining each of the 
terms of eq (1), Blackadar ° [1] pointed out that 

"the independent estimates of Mx and Aly did not 

agree and negative values of each were numerous and 

frequently large. It is assumed that the reason for 

the poor performance is in the inaccuracy of the 

geostrophic wind." 

For this reason, he devised an indirect method for obtaining the 
values of Ug and Vg which would minimize the sum of the absolute 

‘magnitudes of y= My at the four levels (950, 900. 850 and 800 mb). 
He forced Mi, to be equal to My at a single level and at the oe 
standard ee oeeatts times; then the geostrophic wind components were 
computed. Next utilizing the thermal wind relationship and assuming 
the vemperature sounding to be accurate, the values of Ug and Vg were 
found in the vertical, The variation of geostrophic wind at times 


other than the standard radiosonde times was determined on the 
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intervals and three hompiy S48, ..tervals. Tne valiies of rand @atE 


SO obtained were cenerally a tcyog@Ment at all times, 2lackadar Li 


aecided 
Fe eov0 Gispose of the remiining differences, the best 
estimate = was assiried to be given by 
ih 
? 
Fis. 3 shows Blacksdar’s | 1° cemputed values of ™ as a function of 


merrhc and time. 
\ 


It is desired to find a fiaction (‘w (z,t) which will yield these 


fois of MY . Due to the rerid@aie raturg of the solittion, fir. 3 
must be modified in order to show the same values of #1 Teor 


corresponding heights et the berinning and end of the 24 hour period, 


Pre alceration ef -rir, 3 was minimized insotar as possible, and the 
modified fig. 3 ia shown as fig. 41. 


The followin: @pproach Vo the problem was used. Determine the 
Merc series which fita the observed values of @ far i2 tines 
( 2 hour intervals) at each of the given heights 0, 200, 400, 600.... 


1E00 meters. An exemple fur the height cf 200 meters is 


Pro sAa(2,t) = A. (2 ye AY (2) cos ince + Ty (2) | 


me ‘| ic =u (4) 
ea, (2) cos | (2 # T.(2) | #Bcoedg cos [ (6t + Tg (2) 


"ere Sine and cosine terms belonring to tne same harmonic have been 
combined: A, (2) is the mean value, A,(2) is the amplitude of the 
first harmonic, T, (2) is the puasge ancle of the first narmenic in 


decrees end t 15 time in huurg efter +t te Gisele Coli exoressed in 
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gevrees) Vow, having 10 ecuations of the same form as eq (4), fit a 
continuous function of z to the 10 values of A and T for each harmonic 
Woich males a sienificant contribtizion, This would yields ne final 


moak-Ol ne eCdue ~ 1 OMe 


MM (z,t)= Ag(z) + Ay (Zz) cos & + 7 (2) | 


+ Ao(z) cos | ot 2 To(2° | frome) (COs |6¢ + 26(2)| ©) 


The hermonic analysis at each even 200 meter height was accomplished 
yielding the values in Table l. 

Fig. 5 shows values of Ay (the mean), Aj and Ap versus height; 
fig. 6 shows Az, Ay and dg versus height. Fig. 7 shows Ty and T9 
versus height and fig. 8 shows Tz, Ty and Ts VeneuGce Gel pic. 

@ diagonal difference table was computed for Ag values versus 
Mem@entce., ithe second differences of Ay, AA, showed the most nearly 
constant values, Although the second differences were not very close 
(ranging from 45 to 128), a parabola was fitted to the values of Ao 
by the least squares method. As expected, the parabola did not fit the 
data a Wienin Sufficient: accuracy. 

Next it was noted from figs. 5 and 6 that a possible orderly 
Meerease in whe ratios of A. /Ko existed at each height as i increased 
wecme! to 5. A plot of the ratios A, /Ap versus height was drawn 
(see fig. 9). Utilizing these ratios the equation for M (z,t) would 


taize on the forn 


Fl = Ag (2) oa +B, (z) cos [t + T (z)| + Bo(z) cos | (2 + To(z)| 


; (6) 
eee va) .cOs | (6+ + Ts (2\} } 


te 


TABLE 1. 
Results of harmonic analysis at 10 levels. 


Z is height above surface (m), A, the mean coefficient of eddy viscosity 


= O i = 
(gz cm Leecly, A, the amplitude of the i&® age (g cm 1 Sec my and 


ry 2 ¢ ° 
r; is phase angle of i” harmonic degrees. 


Ay ky Az | T3 Ay Ty Ac Ts N6 


40.0) 5407 |-22.9| 754.91) geomeoa.e| 226 
35.6161.2] 24.4] 83.7] 19,.4/125.8] 8.8 
OHS) 69e |) lonol lla. > [elon e560) 9.3 
19.4)69.4} 15.3) 123.8] 2%3.7|165.1/12.2 
Hoecloa7| 18.5) 1065 |) 2oment 52.21 14.4 


1000 ‘ +5] 16,01158.11152.6) 95, 0\eo. orn OGe 18.7) 31.2 
1200 : : 65| 8.21351.1| 8.3/248.4) G6.25119626)) 6.9] 93.7 
1400 : He sere 8.31360 8.0 |258 8.0} 114.6 
1600 8.6]756 111.8] 438 15.6 1266 5. 211105 





1800 | o4 18.1|984 | 31.0] 539 Ah Oo} 249 45.5160 


(1) All surface values are the result of extrapolating Blackadar's [1] 


curves down to the surface, 
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Yer? An(z) 1s tne function of height needed to fit the observed values 
of the mean Mat each height. By (z) is the function of height needed 
Zot wie rae Lo B,/ 40 at each height and Ty(z) is the function of 
height needed to Pit the observed Ty values at each height. If the 
curves of fig. 9 could be fitted by a relatively simple function of 
height vice an estimated 5th or 6th degree polynomial, there existed a 
possibility of a series analysis and possible further simplification. 
One noteworthy aspect of this approach should be observed. + appears 
that a fairly simple function of height (e.g., a parabola) could be 
fitted to the ratio curves (fig. 9) with reasonable accuracy if a 
iayerized approach is followed. The dashed lines on fig. 9 represent 
the probable curves which would fit a 2 layer model. The division of 
layers appears to be in the vicinity of 600 to 700 meters. This 
conceivably could be of physical significance. 

Next it was noted that A, (fig. 5) increases non-linearly with 
height above 460 meters. Further, the point z = 460 meters, for which 
= T3g.om™*sec™*, appears to be a point of symmetry. The axis was 

O 
translated to this point and Ap versus height was plotted on log-log 
paver; this plot coincided almost exactly with a straight line. 


n 
By the least squares method, Ap(z) = az + b was determined. 


Tne equation for do WAS 
A(z) = 6c [2-46 [22 4 75 (7) 


where 2 is height above the ground in hundreds of meters. Table 2 


$9 


shows a comparison of the computed values of Ay and observed values 


(from the harmonic analysis). 
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Sinilarly A, (2) and A,(z) were obtained. A comparison of A, (z) 


computed and observed is shown in Table 5, and for A, (z) in Table 4. 


tt 


The equations developed by the least squares method are as follows: 


opal 


W 
> 








Aj (z) = 68 Z- 4,5 + 81 (8) 


Il 
ye) 
Uy 


A, (z) |z - 4,0 |e ee) (9) 


When locating the approximate point of symmetry of the Ay versus 
height curve (fig. 5), it was necessary to employ the method of least 
Squares to determine which of three logically chosen positions for the 
heisht axis yielded the least difference between A(z) calculated and 
A(z) observed. In spite of this effort, relatively large differences 
were observed in the range, surface to 400 meters (see table 4), where 
the symmetry was weakest; however, the fit appears to e good above 
500 meters. 

A linear approximation was used to fit T,(z) and T,(z) (see ig. rae 


The eguations determined by the least squares method are 


~ 1.20: 2a ee 


Ty (2) 


(10) 


LO) ~ 1.81 2 + 39 


where z = height above surface in hundreds of geters and T, (2) and T, (2) 
are in degrees. Table 5 shows a comparison of the phase angles computed 
by eq (10) for the first and second harmonics and the observed phase 


angles. 


-_ - 
- — 
‘ =—= © “Sea bo@ 
- - - -— - 
ae 6 _ 2 a 
— —— —_ @«& 
-— == & = 





TABLE 3. 


A, (z) (g sm sec7!) computed versus A,(z) observed 


HEIGHT COMPUTED ODSERVED PERCENT 
Le 0 Ay (Z) A, (2) ERROR 
0 25.9 

2 ey 

th TS 

6 oe 

8 5.8 

10 8.2 
12 4.3 
14 vey 
16 463 
18 fal 
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TABLE 


ho(z) (g Pie aeocn) computed versus Ao(z) observed 


HEIGHT COMPUTED OBSERVED PERCENT 
reer seh: (a) A, (2) ERROR 
0 110 86. Pa SS 
2 48, 74. 24, 
4 D2 Dhar BOD 
6 48, 50k 270 
8 110. 98. Liles 
10 2s 20D. 2.4 
2 450. 471. 4.8 
14 one 806 6. 
16 1157 1169 Lee 
18 1674 1620 | oa. 
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TABI. 
T, (2) computed versus T,(z) observed 
and To(z) computed versus To(z} observed for 10 levels. 


T, (2) and T(z) in decrees. 








T, (z)obs. Ty (z)comp. To‘z)obs. To(z)comp. 


O 20%: eae 40, 39 
c ee 18. 25 35. 
dy 13 16S Bek 72 


ea 








Jnon examination of the remaining amplitude coefficients and their 


corresponding phase anrles (see figs. 6 and 8), it appeared that the 


+ 


t 


LO 


functions of height needed fit these values to the desired degree of 
aeecuracy would be relatively long and too complex to warrant further 
watore. an thes piace. 

Also one of the objectives was to see how complicated a function of 


heicht and time it would take to yield the observed values Of Aas elects Ds 


whe equation for MM (z,t) DAMS ear SOU ed, 1S 











M (23%) = 6.0 gn A 6 or) 4 73 4+ (8. 4 Z- 4.5 e + 61| 
cos “T= (ee 2) |e [2-401 z ESS 55 (11) 
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after t = 0 (1500 CST). Values of M (z,%) were computed using the 


Ww 


Marchant desk calculator for time intervals of one hour and every 100 
meters heisht. Fir. 10 shows computed values of M plotted on a 
Meieht-time coordinate sysicn. 

Comparison of fig. 10 and the desirod M distribution from fig. 4, 
snows fig. 10 to contain a fairly representative set of isolines for M 
from 1000 CST to 2000 CST from the surface to 1800 meters. Both 
Figures show the largest order of magmitude for the coefficient of 
woe Viscosity at approximately 1500 Cs? et all Vevels. In addition, 
bovh figures show a minimum for {in the vertical at approximately 
400 meters between 1000 CST and 2000 CST with M decreasing from the 


surface to 400 meters throughout the period and with Af increasing 








"SDIUaM AMY OM]. ysarp 
am pu» uvaly papysd Puy. hg ‘4 ae) 4ybray PYyY uw do uo oun? 
Se uoprquysip hysorsia hppa Jo 4uarda02 ey) 'Q]'S| 4 


Susi Dbo JEO D12 Lie 
Z2bO Z £0 Z\¢ 


TS Nie wee? 
mo NZ UZ VN Zen 
ya SA ss I NE Zak 

—~ WWE TIZ SAI 
om \ We NUE MLE Zs 





2 


Oo 





: a : | 909) 
ot \ MT NULL LA \\WWe WHEL Ze 
= ( _ —- a os x FT 008] 
vi xv | Cee & 7 Re OQ FO 0 ~ Ay 


(ssajau) |LHOISH 


24 








Ee | aA 


reavidly with height above 400 meters from 1200 50ST to 1900 CSTs In 


Me 
af 


co 
a 


. 10, however, the decrease of M with time at any level above 1000 
Revers 18 too great in vhe vicinity of 1900 CSi es tome ceinercasce of 
MM with time too large at any level above 1000 meters in the vicinity 
of 1000 CST. Also, the vertical fradient of M above 400 meters at 
1500 OST is less than that observed on fig. 4+. The latter could be 
rectified by computing the remaining harmonics which would tend to 
increase the values of A computed at 1500 CST (t = 0) so that they 

more closely approach corresponding values of observed M (Cigale). 

Between 2000 CST and 1000 CST, surface to 1800 meters, a 
comparison of fir. 4+ and fig. 10 shows fig. 10 to have a completely 
unacceptable representation of the distribution of M, Pig. 10 snows 
aecolumn of minimum M values extending vertically from the surface to 
1800 meters at about 2200 CST, with no phase lag, and another minimun 
column, extending from the surface, at 0500 CST, to 1800 meters, where 
the minimum occurs at about 0700 CST. These minima are separated by 
a relatively strong maximum, surface to 1800 meters, occurring at 
0200 CST. In contrast, the desired form, fig. 4, shows a strong 
single minimum extendins from the surface, where it is estimated to 
occur at 230C CST, to 1800 meters where it is observed at approximately 
O80 CST. 

The double maxima, double minima vice the desired single maximum, 
Single minimum, are attributed primarily to the fact that the mean and 
only two of the six harmonics were utilized in deriving the values of 
M for fig. 10. Fig. 11 is introduced to lend credence to the above 


remarlx, This figure shows the resultant of the first two harmonics 


versus time for height equal to 1200 meters, The resultant M 
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disuribution utalizine tie mean and the Tirstl Guo Pmarmom@sca as 
determined by the harmonic analysis, not from fitted values, tends to 
show the following characteristics3an absolute maximum at 1500 Col, a 
relative maximum at 02350 CST, and two minima, one at 0/00 Cat and the 
Gener on. 2200 Cat, 


Bach of the above maxima and minima are observed on fig. 10 at 


the corresponding times and the magnitude comparison is as follows: 


TIME fo erojaemogi yes ju) JM fron gap oe rae 
1500 CS? 1630 1620 
2900- CST ~85 Sil 6, 
©2350 Csr ALS oho 
OF 00-CST a iew ~140 


Similar analysis at other heights would yield similar results. 

This would tend to show then that the fitted values have not 
distorted the picture of M distribution with time and height except to 
a slight degree and that the cause of the deviation of the computed M 
values (fic. 10) from Bleckadar's l1 | observed values (altered 
fic. 4) must be the remaining four harmonics. Further, when using 
eq (11) and including the remaining four harmonics, the values of 
JM computed compare very favorably with the observed values. 

Feruinent to the variation of M With height, S3stoque a 
has recently approximated the distribution of M by a linear decrease 
of M above 4 maximum value found at a height of about 50 meters. 

i@ proposes a two layer model for the frictional layer, the lower 
layer, called the surface sub-layer, extends to approximately 50 
meters, and an upper layer, called the transition zone, extends to 
epproximeately 2000 meters, Utilizing this §{ distribution, he has 


numerically integrated the ecuations of motion and heet conduction 


et 
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equation to obtain the 24 Nour distribution of wind and temperature. 
ime Czurnel vind €istricuvion hus ovvained O77 uacoduc [ 3 appears 
~O- De IN Palrly- food arreeméent. with che diurnels gimdecharsc eri scics 
observed on several occasions. However, for Elackadar's [2] dace. 


}, it was necessary to attain an averare height of 


AN 


(see fi 


ys 
Coa 
approximately 1150 meters before A/ varied linearly (but increased) 

with height to within a 10 percent error of the observed value of H . 


Purcher, the slopes of the lines representing the linear variation of 


MM with height (for different times) varied up to 55°. 
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ae 


fhe Coefficient of sdde “Biscosity as @ Function of Stanilivy and 
onear of the wind. 


dind date Lor this phase were obvaised -fromerw7sc. ) 250, 
reproduced from Dlachkadar Ea - Temperature data for the 24 selected 
days were obtained from photo copies of WEAN~43 for the months of June, 
Jul:,, Ausust, and September 1952 for Vichita, Kansas and Oklahoma City, 
Olclahoma, which were furnished by the National Weather Records Center, 
Ashville; North Carolina. 

All tenperatures at the revorted levels (surface, 950, 900, 850, 
800 and 750 mb) at both cities and at both sounding times (1500 GOT 
and 0200 GCT) were converted to potential temperature. The surface 
potential temperature was obtained by the Arowargram and upper level 
potential temperature was obtained from the Smithsonian Meteorological 
Tables. 

For the 0300 CCT sounding and at each level, the average potential 
temperature was obtained by adding the 24 observations at Wichita to 
the 24 observations at Oklahoma City and dividing by 48. This yielded 


a 


the averace sounding for the two cities over the 24 day period at 0300 
GOT. <A similar procedure was followed to obtain the average height of 
the 950, 900, 850, 600 and 750 mb levels. 4All heights are in terms of 
0.95 ceodynamic meters, The same procedure was followed at 1500 GCT. 


the averace heights and potential temperatures obtained for the 


two times are showm in table 6. Next the stability factor 


wes computed for the center of the averaged layers using finite 


erirerences [orszhe parcia. derivavive of O . ovebllicy factor 





TAB GeO. 


Twenty four day average heichts (in 0.98 reodynamic meters) 
and potential temperatures (degrees K) for the reported 


levels of concern. 





0300 SCT 1500 GOT 

level . Te © 
degrees K meters derrees K 
Surface % Surface 0 HONS 
o50mb 505.0 Rinks 201 502e 
900mb 506m 900mb Gere | es. | 
850mb 307. 850mb 1168.2 50 fee 
800mb KOSe Ss: &00mb LoSinO 308.6 
75D0mb 


311.4 750mb 2237.6 | 310.7 











versus heicht of center of layer for the two times 2s shown in table /. 
; aes 
du ov 
Values of the components of the vertical wind ee a ) 


for each 100 meter level were then obtained from fics. 1 and 2 


utilizing finite differences for the partial derivatives. The shear 
7 Ou LL ww \* | G 
ractor, | (24) {x [O was computed next for each 100 meter 
oz Oz / 
evel tor the wilco times. 

Graphs of 4 observed versus the stability factor and of M observed 
versus the shear factor were then constructed in an attempt to find a 
correlation between these factors and M « One apparent correlation was 
that in the lower eee (below 600 meters) where pronounced stability 
variations occurred, strong stability was associated with a small Ve , 
end vice versa, Above about 600 meters the contrast in stability factors 
sasmconsideréd neglifible. Moreover, there wes insufficient difference 
between the M values at each level for these two times (0900 CST and 
2100 CST), nor between the shear factors in the upper levels (above 
600 meters), to detect any correlation between wind shear and the 
values of | observed. In the lower layers, when the snear appeared to 
be larre, the value of was also large. In order n@ OE ss tae 3 
correlation between wind shear and | existed in thesupper levels, it 
was necessary to examine the data at times of minimum and maximum MH 
(i.e. 0300 CST and 1500 CST). Computation of the required shear factors 
shoved a strong correlation between low shear values and high 4 values 
ome vice versa for all levels. This did not contradict vrevious resul¢s 
above 600 meters but it was in opposition to the relationship found 
between shear and below 500 meters at 0900 and 2100 CST. 

Since an expression throughout the vertical for M in terms of 


% ° 


Sstapility and wind shear is a most desirable form, a compromise eppecred 


pal 





TABLE 7. 


Beloit, toactor, =. 
or 


24 day average soundings at 0300 GCT and 1500 COT. 


»versus neight of center of layer for 


| oe, 1" 
z 


0300 CCT 


HeLent. of 





Layer center af Layer 
Lawexr 
\meters ) 

Surface to 950mbk 94, 7498 [Surface to 950mb 100 0807 
950mb to 900m 424 0826 950mb to 900nb 447 1949 
900mb to 850m 908 0865 900mb to 850mb 920 1265 
850mb to 800my 1417 - 1065 850mb to 800mbl 1428 0969 
800mb to 750m 1953 Pio 800nb to 750mbh 1963 es 6 
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to be necessary, i.e., a two layer model with Ma function of the 
stability factor below 6500 meters and M a fincelon ef the snear faccor 
above 600 meters. Now, from visual inspection of figs. 1 and 2 at 
1500 CST and from prior shear computations the shear can be seen to 
be very small and practically constant throughout the layer surface 
to 1800 meters; yet the value of /4 (see fig. 4) increases rapidly 
from 600 meters upwards. Therefore, a second compromise from the 
originally intended solution appeared to be in order, i.e., above 
600 meters let M vary with the shear factor and with height. 
Further, due to the lack of temperature data below 600 meters, it was 
decided to utilize the mean M (z) with respect to time in an attempt 
to more accurately describe the M observed values and also perhaps to 
extend the coverage below 600 meters to the entire 24 hour period 
rather than merely the times of the two temperature soundings. 

The general approach used in the solution of the problem is as 
follows: 


1) Assume (for both layers) 
A = Fay t 44 (a2) 


where H is the observed value of M for any point in time and space, 
M (z) is the mean value of all 1 observed at the level z, and Ay 
is the difference between the observed value of Mf and the mean value 
of Mat that height. 

2) Find the best function of height which yields the A (z) 
aestri bution. 


3) Account for AM , the deviation from the mean, 


WN 
\N 





a) below 600 meters by a function of stability factor, 
b) above 500 meters by a function of shear factor. 
) ed 


A function was found in Section 2 to represent the mean A{ at any 
height (see eq (7)); however,this was based upon values obtained from \ 
the isolines of A{ in fig. 4. For this section an attempt was being 
made to find the function which yielded the computed values plotted 
on fig. 3, inasmuch as they were in part computed from the vertical 
shear of the horizontal wind. A similar approach to that used in 


deriving eq (7) was in order and yielded the following for (4 (z) 
Fi (2) = 5.6 |2 -4.8 laa a (13) 


Where Z is height above the ground in hundreds of meters. Due to the 
Belacively ¢ood fit, this equation was taken to apply in both layers. 

hex: AM values were computed by subtracting the M (z) values 
computed with eq (13) from the values of Mobserved on fig. 3. 

A plot of AY versus spabi lity (acvor dseanevmeaset1¢. lee 
Four of the points plotted are from observed data for the period 
(table 7) and the fifth point is an estimated value for a height of 
150 meters eve, ILS O16 16 sye This value was derived from 12 soundings 
@eeen during the Great Plains Turbulence Field Program bi in 
August and September 195% and under similar atmospheric flow 
conditions. The stability factor obtained was .0035 ae. and the 
AM corresponding to 1500 CST and 150 meters is Hog cm7tsec7, 
The curve appearing on the figure has been faired in and does not 
represent the Pitted function given below. 

Cr various functions of the seabilivy [ecler- tried. ciat 
which yielded the smallest "least squares" error for Aa was the 


hyperbola, 
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| is 
-121 = 5 (S—) +% |200,560 ~ 7400(39) + 127 (39)? |* 
< see 
where (3, ) Ei ois X io (14) 
| 0 0 


Combinine eq (13) and eq (14) then yields the desired form for 


Senvucine (4 5 ise. 


J computed = 5.60 | Zz ~ 4.8| 2-20 + 186 -5.€ (Sg) 


(15) 
i 
+a [200, £ 90 ~ 7400(Sg) + 123 (Sq)2 | 


wa m2 
Here Zz is in meters x 10 -, 


Table & shows a comparison of M computed and / observed for the 
layer below 600 meters, The largest error occurs at 2100 CST and a 
height of 424 meters, It should be noted (see figs. 1, 2) that at this 
point the shear is relatively strong which, according to eq (3), would 
mena to yield a small value for AY. aA similar suaieiemcites for the 
next larrest error, percentwise, at 0900 CST and 437 meters. Here, 
where the.observed value is greater than the computed value, a low 
value of wind shear is observed. From the preceding analysis, it 
apvears that shear is becoming a factor in the upper levels of the 
lower layer. In this transition zone from approximately 400 to 600 
meters, a better representation of M might be obtained by considering 
shear as well as stability. 

For the layer above 600 meters, a graph was drawn showing height 
versus shear factor and the associated A value (next to each point) 
needed to yield the observed value of M for the four times shear 


computations were made. This is shown as fig. 13. 





TABLE 8 ° 


a 


> 2 =) =] 
Comparison of Af computed and otbserved(g om “sec 
x : b] 


surface to 600 meters. 


Time Height above Surface /M observed FA computed. 
ccor (meters) & cm -sec™ Z@ ca ~sec™ 





150 hso 450 
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AS “reViousiy Svc CBR. jas desired to Sccounus or Espa valucs as 


a. function of height anc shear. The approach to the problem was 


Similar 2o that used in section 2. An cquilateral hyperbola: was tapece 
to the observed values of shear and A h level (S00, 800, 1000 
BO ume OuServed yeiues of shear and: JAM av-eac eve O00; Po 


evecol800 meters) to provide the AP distribution as a function of 


Po 


Shear factor. The ré@sult is 


= 
mit : Ey “ 06) 


(SS 
where ie is the difference between the observed M and M (z) at 


height z2 above the ground; c, and bz are the constants determined by 


the least squares method at the z level. Fig. 14 shows c, and b, 
versus height. The curves appearing on the Pigure have been faired 
imoeeanc do not represent time fitted function given below. 

next the c, and by values determined at each height were fitted 


fem COncinuous funccion of height. The functions of Herehl ee 


revresentins each of these is 


c(z) 


il 
e 
ae) 
MN 
o™~ 
N 
~~ 


(17) 


We ani. a (18) 


~2 
where Zz is heizht above fround in meters x 10 ~, 


A comparison of the c, and b, values observed and the c(z) and b(z) 
velWes computed is shown in table 9. The resulting function of shear 


and heieht representinz the Ak values at 600 meters and above is 


AM 22 @  - 6) (27> (19) 
yU\> Lf av 
Gay Sy f° 
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The constants c, and b, versus c(z) and b(z) at each level 


e(z) 





A. 7457 9381 
16 14,372 


18 of 147 


a 








and the function representing M in the upper layer is 


XN 


MM 25.6 (2 - 4.8)2% 65 + 22 yo Ci) 
) ov ) fee 
oz} (doz (20) 
A comparison of the computed M values and the observed values is 
shown in Table 10. 

The distribution of M computed indicated in Table 10 compares 
fairly well with the observed values. The salient features of the 
observed H distribution are preserved by the computed pattern, i.eé.; 

a maximum at each level at time 1500 CST (one exception at 600 meters) 
and a minimum at each level at time 0300 CST. In most instances the 
order of marnitude of the observed and computed values are the same 
(five exceptions out of 28 observations). The value observed at 

0900 CST, height 600 meters, is approximately one fifth that computed 
by shear and height function. This discrepancy may partially be 
accounted for by the higher stability factor observed in the vicinity 
(see table 7). The negative values computed for heights 1200 and 1400 


meters at 0300 CST are primarily a result of the curve fitting process. 


(See fig. 14 and table 9). 
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TABLE 10 
Comparison of computed coefficient of eddy viscosity, /((g on” + see7!) 


and observed for four times in uvper layer, 
+] Oo uv 


Height 


(neters) Mc Po ) Me es) Jc Mo) Vase (Ho) 


| 1800 | 7607 |(7020) 455 |(480) 105 (328) 
| 6709 |(5900) 7344352) 36 {(97) 309] (208) 
| 5800 |(4680) 2200} | ~29 |(65) 
| 1200 | reir |(30s0) Mais) =7 (37) 120] (112) 
| 1000 | a6 ftr250) 4 |63) 7 (81) 
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4. Conelusions,. 

From the approach followed in this investigation a simple 
representation of the diurnal variation of the coefficient of eddy 
Viscosity with heignt was not apparent. 

I% has been show that the function of height and time 
representing the mean and first two harmonics is inadequate to show the 
pulrene Peavuures 101 —vme observed distribution, 

In order to represent the observed distribution of the coefficient 
of eddy viscosity by a function of height and time, it would take a 
relatively long series of terms, estimated to encompass at least five 
harmonies, involving exponential functions of height, fourth or fifth 
depree polynomial functions of height, and trigonometric functions of 
angles Pe eient on time and funetions of height. 

A two layer model is developed in an attempt to account for the 
@eserved cistribution of the coefficient of eddy viscosity as 2 
function of stability and shear of the wind. It is found useful in 
che lower layer, and essential in the upper layer, to include height 
above the ground as one of the parameters determining the distribution. 
The variance of the coefficient with stability in the lower layer 
appears well founded, as does that between shear and the eddy 
Viscosity coefficient in the upper layer; consequently, a fair 
representation of the general characteristics of the observed diurnal 
variation of the coefficient of eddy viscosity is obtained from the 
empirical formulae. representing the two layers. 

Purther imamestigetion uvilizing more frectient and detailed 
and temperature observations would be likely to yield better results 

; 


and might eliminate the need for the two layer approach. 
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